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We benchmark three standard approximations for the many-body problem - the Hartree-Fock, 
projected Hartree-Fock, and random phase approximations - against full numerical configuration- 
interaction calculations of the electronic structure of atoms, from Li through to Ne. These 
configuration-interaction calculations used up to 2 x 10 s uncoupled basis states, equivalent to 10 7 
coupled basis states (configuration state functions.) Each method uses exactly the same input, i.e., 
the same single-particle basis and Coulomb matrix elements, so any differences are strictly due to 
the approximation itself. Although it consistently overestimates the ground state binding energy, 
the random phase approximation has the smallest overall errors; furthermore, we suggest it may be 
useful as a method for efficient optimization of single-particle basis functions. 



I. INTRODUCTION 



In principle we know how to find the electronic structure and binding energies of atom: solve the many-body 
time-independent Schrodinger equation for N electrons about a nucleus with charge Z: 

^= \lt-f^- — +EAr = M w 

\pt 2m n ^ r, r,, J 

or variants thereof including relativistic corrections. In practice one cannot solve this exactly, or even fully numerically 
much beyond helium. 

So one turns to approximations. But for any approximation there are errors intrinsic to the approximation method 
itself, and errors due to the input, such as truncations in the space or assumptions in the Hamiltonian. Surprisingly, 
rather little effort has been put into disentangling these two sources of error. 

In this paper we benchmark three approximations against an 'exact' calculation with a fixed space and fixed 
nonrelativistic Hamiltonian. For our exact calculation we use full configuration-interaction ( full CI); for a chosen 
single-particle basis, we can find numerically converged solutions for the ground state (as well as some number of 
excited states). The full CI spaces are large, up to dimension 10 s in the uncoupled ('M-scheme') basis, or 2 x 10 6 in 
the coupled basis of configuration state functions (CSFs); see below for further discussion. We also have a suite of 
codes that carry out Hartree-Fock, angular-momentum projected Hartree-Fock, and random phase approximations, 
using the identical input as CI. Thus our goal is to carefully document errors due the approximations themselves (and 
not due to the single-particle space, relativistic corrections or lack thereof, etc.); in this we are extending prior work 
benchmarking approximations for the nuclear many-body problem [J-Q- 

In the next section we define our many-body methods and discuss briefly computational issues. In Section 3 we 
give our results for several different spaces for second row atoms, from Li to Ne, comparing not only the ground state 
binding energies but also the experimentally more relevant ionization potentials and electron affinities. Because full 
CI is very time-consuming and our approximations computationally cheap, in Section 4 we consider the possibility of 
using approximations for fast optimization of the the single-particle basis, and of the considered methods find RPA 
the most accurate. 



II. MANY-BODY METHODS AND INPUTS 



We work entirely in occupation or configuration space, meaning we start with some set of orthonormal single-particle 
states, {4> a (f)}, with good angular momentum and parity, and introduce in second quantization 0,0 fermion creation 
and annihilation operators a\a that create or remove a particle from the states <j>(f). For definitiveness we have N s 
single-particle states. The Hamiltonian is then [IHH 

H = H lb^b + \Y1 V ab ,cda\a\a d a c (2) 

ab abed 

The one-body and two-body matrix elements, H* b (which includes both kinetic energy and the nuclear central po- 
tential) and Vab.cd respectively, are appropriate integrals over the Hamiltonian and the single-particle states. Because 
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the Hamiltonian is an angular momentum scalar, we generate the matrix elements in coupled form. Given the single- 
particle space, the matrix elements are computed externally to all our codes and read in via a file. Any radial form 
of the single-particle states and any form of the interaction, including non-local interactions, are allowed and make 
no practical difference to our many-body codes. Because we work in second quantization, the matrix elements are 
antisymmeterized and we do not separate direct from exchange terms. 

For this paper we restrict the Hamiltonian to the simple nonrelativistic Coulomb interaction (JIJ. We build our 
single-particle radial wavefunctions from Slater-type orbitals (STOs), that is, functions of the form r n exp(— fir), which 
allows us to compare to other work as well as making the integrals analytic. We orthonormalizc our single-particle 
basis so that the one-body part of ([2]) is diagonal. For our specific choices of parameters for our basis, we use the 
CVBn bases sets Q, described below. Finally, we describe the atom as a two-species systems, with spin- up and 
spin-down electron occupying distinct orbits. We fix M s = (n-j- — n^)/2 in all calculations, as we have no spin-orbit 
coupling. 

These methods were chosen because we have existing codes for them, and thus leave for future work other methods 
such as the multi-configurational Hartree-Fock approximation [1, 0] and many-body perturbation theory [3, H, 



A. Configuration-interaction calculations 

For full configuration- interaction (full CI) calculations 0, [H, , we use the BIGSTICK code [17| , developed 

originally for nuclear CI calculations but which can be applied as well to atomic calculations. The only assumption 
is that the single-particle basis be states of good angular momentum and parity. 

BIGSTICK creates a many-body basis of Slater determinants {|a}}, constructed from the single-particle basis: 

N 

\a)=l[al\0), (3) 

i=l 

where N is the number of particles (occupied states). Furthermore, because we have no spin-orbit coupling (a 
generalization we leave to future work), we fix the number of spin- up and spin-down electrons and thus fix M s = 
(n-j- — nj.)/2. Because single-particle states each have good l z , it is almost trivial to generate many-body basis states 
with fixed total M = L z ; in nuclear CI calculations this is called an M-scheme basis. And because the Hamiltonian 
is an angular momentum scalar, the final eigenstates thus will automatically have good total L (and S). Although we 
have the capability for various truncations on the many-body basis, such as allowing only single- or double-electron 
excitations, we make no such truncations here as they have no correspondence for our approximation methods. 

Most atomic CI programs such as CIV3 [l5[ and GRASP92 [l|| use a basis of configuration state functions (CSFs) 
coupled up to good L and S (or total J, hence in nuclear structure physics this is called a J-schemc basis). Such 
a basis is smaller in dimension, by a factor of 10 to 100, but there are always trade-offs: the basis states are linear 
combinations of M-scheme Slater determinants and hence the many-body Hamiltonian matrix elements are more 
complicated to calculate. And although the dimensions of an uncoupled M-scheme basis are larger, the Hamiltonian 
matrix is much sparser. Furthermore the simplicity of the uncoupled basis allows us to avoid storing the many-body 
matrix element, and instead recompute them efficiently on the fly [T^, which further reduces the storage of the 
Hamiltonian by a factor of 10 to 100. We can go up to uncoupled basis dimensions of 10 8 on a desktop machine (and 
similar programs on parallel computers reach up to 10 10 ), as listed in Table II, while the dimension of the largest 
equivalent coupled CSF ground state dimension is about 10 7 , as listed in Table III. 

Given these arguments, it is interesting to compare basis sizes. These are not always given in the literature, 
and many 'large-scale' calculations rely upon careful pruning of the configuration basis [19h21| ]; nevertheless, for 
comparison purposes a lar ge s cale relativistic CI calculation [22( used up to 197,426 basis CSFs, while a more recent 
nonrelativistic calculation 23[ used up to 901,816 CSFs. Quantum chemists surpassed the one-billion determinant 
limit some time ago [3 [H ] . 

In a given uncoupled basis, BIGSTICK computes the many-body Hamiltonian matrix elements (a\H\fi) and finds 
the low-lying eigenstates, including the ground state, using the Lanczos algorithm [25j (the Davidson-Liu algorithm 
[2IH [27] | is efficient for the diagonal-dominated atomic problem but not the more general many-body problem such 
as nuclear CI; furthermore recomputing the matrix elements on the fly is better suited for the Lanczos algorithm). 
Convergence of, say, five states take a few hundred Lanczos iterations; on a desktop machine with eight OpcnMP 
threads we can accomplish this for an M-scheme basis in a few hours or at most a few days. 

With the single-particle space and Hamiltonian matrix elements fixed, the resulting eigenenergies for the full-space 
CI calculation are numerically exact. All other calculations described in this paper are numerical approximations to 
these results. 
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CI calculations have been previously carried out for second row elements, see for example nonrclativistic [2814341 ] 
and relativistic calculations [22l[35j , though most of these were not full CI and thus were in much smaller dimensions. 

B. Hartree-Fock approximation 

The Hartree-Fock approximation (HF) is a variational method using a single Slater determinant iMll. In 
occupation space @ one applies a unitary transformation among the single-particle states: 



N 

a=l 



\=Y,U ia al (4) 



and creates the trial Slater determinant 



N 



i*T)=n e ii°)- ( 5 ) 



Then one finds a Slater determinant that minimizes the energy, that is, that minimizes 

E= (lll^M ( 6 ) 

Because our input matrix elements are antisymmeterized, we fully include both direct and exchange terms. Working 
in occupation space the exchange term causes us no difficulty (or, to put it another way, any difficulty is off-shored 
into calculation of the antisymmeterized integrals). 

Many Hartree-Fock calculations enforce good angular momentum, for closed-shell systems or closed-shell plus or 
minus one particle. Our HF code allows the Slater determinant to break rotational invariance, even for closed-shell 
systems; such Slater determinants are 'deformed.' Experience in nuclear physics suggest this to be a fruitful path. In 
a few cases, when possible, we will compare and contrast spherical and deformed HF solutions. 

In order to minimize, we use the standard Hartree-Fock equations, which consist of iteratively solving 

hui = erfi, (7) 
where the Hartree-Fock effective one-body Hamiltonian is 

hab = T a b + 2^ Vac,bdPdc (8) 
cd 

and 

N 



Pdc 



J2 UicUid- (9) 



Internally we use an N s x TV rectangular matrix VP, with the N columns representing the transformed occupied states 
cj. In that case p = St"* 1 . 

(Because we leave out spin-orbit coupling and conserve n-j- and nj., we use separate Slater determinants for spin- up 
and spin-down electons; this is a straightforward generalization. ) 
The Hartree-Fock energy is now 



N ( 1 

( £l ~ 2 ^ Vic ' id 

i=l \ cd 



Ekf = > c- - 77 V V^dPdc , (10) 



where we have to subtract off the potential to keep from double-counting in the single-particle energies. 
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C. Projected Hartree-Fock approximation 

If one does not have a closed shell, a single Slater determinant will not in general have good angular momentum. 
Conversely, an open-shell state with good angular momentum, which can be written as a sum of Slater determinants, 
will have important correlations that lower the energy. 

One can start with single-particle states which have good angular momentum and then construct many-body states 
with good total angular momentum; these are called configuration state functions (CSFs), and have been applied to 
second-row elements, e.g. (37l. l38|. 

We take a different approach, however. Our Hartree-Fock code allows for arbitrary deformation: the single-particle 
states generated by the self-consistent field do not have good angular momentum. The Slater determinant built from 
these states is called a 'deformed' state. In order to project out good angular momentum, we introduce the standard 
projection operator [f| HI] 

p j mw = J dnv$*,(n)R(n), (ii) 

where Q = the Euler angles a, 0, 7 and d£l — dasm.j3dj3d^, P]^, is the Wigner D-matrix [4(|, and R(Q) is the 
rotation operator. To compute {^\R(Q,)\ip) and (^f\HR(ft)\ijj), we use the matrix representation 

i?|*)^R* = * (12) 

where St and arc N s x N matrices representing our original and rotated Slater determinants, respectively, and R 
is a square N s x N s matrix; in our single-particle basis the matrix elements are straightforward: 

Rub = (j a m a \R(a, /3,j)\j b m m ) = $ia.h T) m a ,m h { a iPil)i ( 13 ) 

where TP , is the Wigner D-matrix [|(|. It is useful to note that the j,m,m' for the rotational matrix R are those 
of the single-particle space, while the J, M, M' for the projection operator (fTTj) are those of the many-body space. 
Matrix elements between two oblique Slater determinants is relatively straightforward fill ]. First, 

(^\R\y) = (*|*) = det tftf. (14) 
Calculation of the Hamiltonian (and other) expectation value requires the density matrix, 

, 

Pab (*]*) 

Then 

= H abPab + 5 E V "bMPacP' bd ~ p'adPbc) (16) 
ab abed 

Accounting for two species is straightforward. Note that even if our original Slater determinants were real, by rotating 
over all angle they can become complex. 

We only project out good orbital angular momentum L. We could in principle project out exact spin S; we leave 
such a modification for future work. Because we leave off spin-orbit coupling, however, we automatically get states 
with good S already. 

We now can compute the Hamiltonian and "norm" matrix elements, 

h J M>MI = {^\HP J M ^), n J MtMI = {^\P J M ^), (17) 
which are both Hcrmitian. Then wc solve the generalized eigenvalue equation 

/] ^M,M'9M' = EpHF E n M,M'9M' ■ (18) 
AT M' 

where we allow qm to be complex. Although the deformed Hartree-Fock state will have an orientation, the final result 
will be independent of orientation; we confirmed this by arbitrarily rotating our HF state. 

Both the Hartree-Fock and the projected Hartree-Fock energies will be variational upper bounds to the exact CI 
energy, with the projected Hartree-Fock energy lower than (or equal to, in the HF state has good rotational symmetry) 
Erf- 



I) (i 



(15) 
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D. Random phase approximation 



Once one has a Hartree-Fock calculation, it is straightforward to go on to the random phase approximation, using 
the matrix formulation [l|, Q : 



A 

B* 



B 

-A* 



= Ml) 



Xx 
Yx 



(19) 



where 



and 



c)c n 



H 



(20) 



D 



c]c n ,H 



(21) 



Here i,j label occupied single-particle levels in the Hartree-Fock state, while m,n label unoccupied levels. 

From this we can compute the RPA correlation energy, which is a correction to the HF energy. When working 
with a rotationally invariant Hamiltonian, as we do here, it is possible to have a HF state which breaks rotational 
symmetry; this means one could rotate the HF state in any direction and not change the HF energy. In RPA this 
shows up as a zero-energy mode, which must be dealt with carefully when computing the binding energy [HQ. The 
final result is 



E RP A = E HF - ^Tr(A) + \ Mh- 



(22) 



The RPA energy, while lower than the HF energy, is not variational. In previous calculations for nuclei in valence 
spaces [l[ RPA both under- and overestimated the binding energy; in our chosen atomic cases below, RPA consistently 
overestimates the binding energy. Because we start from a deformed HF state, our RPA calculations do not have good 
angular momentum. (RPA is often said to 'restore' rotational symmetry, but this is not restoration of good angular 
momentum quantum numbers. Rather, the Hartree-Fock energy is invariant under rotation, and this invariance shows 
up as a ft = mode that appears in RPA but not in the simpler Tamm-Dancoff approximation or TDA . In this 
sense RPA 'restores' rotational invariance that was broken in TDA.) 



E. Input spaces and interaction 



Our codes can use (up to memory limitations, of course) any finite set of single-particle states which have good 
angular momentum. For convenience we use the CVB1, CVB2, and CVB3 Slater-type orbitals (STOs) and parameter 
values of Ema et al. Q. These are radial wavefunctions of the form r™exp(— /3r), with the set of STOs and the 
values of /3 optimized for Li through Ne. The single-particle spaces ranged from 8 orbits (for CVBl-Li) up to 22 (for 
CVB3-Ne). 



TABLE I: Basis Set Composition 







Li and Be 






B to Ne 




Orbit 


CVB1 CVB2 CVB3 


CVB1 


CVB2 


CVB3 


s 


6 


7 


8 


6 


7 


8 


P 


2 


3 


4 


4 


5 


6 


d 





2 


3 


1 


3 


4 


f 








2 





1 


3 


g 

















1 


Total 


8 


12 


17 


11 


16 


22 



To construct our many-electron basis states, we fixed rif spin-up electrons and spin-down electrons, with n-s = ni 
or = nj, + 1, which allowed us to regain states of all possible total S. We could and did choose + ^ Z in order 
to compute ionization potentials and electron affinities. We also fixed total AI = L z = 0, which allowed us to regain 
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states of all possible total L. We otherwise had no constraints, e.g., we did not restrict ourselves to single, double, or 
otherwise excitations; although our CI code can make such restrictions, the HF, PHF, and RPA cannot. Dimensions 
of the uncoupled many-body CI spaces we used are in Table HJ while the equivalent dimensions for coupled (CSF) 
bases for the ground states are in Table Hill where we indicated the ground state quantum numbers by 2S+1 (L) 7r . 

TABLE II: Dimensions of the uncoupled (M-scheme) CI spaces 



Atom 




Basis Set 






CVB1 


CVB2 


CVB3 


Li 


192 


1067 


4958 


Be 


866 


10983 


100842 


B 


46116 


909534 


1.1 x 10 7 


C 


300309 


1.1 x 10 7 


2.5 x 10 s 


N 


1.4 x 10 6 


1.1 x 10 s 







6.6 x 10 6 






F 


2.4 x 10 7 






Ne 


8.7 x 10 7 







TABLE III: Size of the equivalent coupled (CSF or J-scheme) CI spaces 



Atom 


2S +\Ly 




Basis Set 






(g-s.) 


CVB1 


CVB2 


CVB3 


Li 


2 S+ 


94 


223 


516 


Be 


1 s + 


190 


775 


3102 


B 


2p- 


7752 


82429 


618245 


C 


3p+ 


35126 


770145 


1.1 x 10 7 


N 


4 s- 


53710 


2.1 x 10 6 




O 


3p+ 


604146 






F 


2p- 


1.8 x 10 6 






Ne 




1.3 x 10 6 







To construct our interaction file, we first orthornormalized the STO basis. We did this by diagonalizing the one- 
body part of the Hamiltonian (and so this depended on the Z charge of the nucleus); then the one-body Hamiltonian 
is diagonal and expressed in terms of single-particle energies. The two-body matrix elements were first computed 
with the STOs, which can be done semi-analytically, and then transformed to the orthonormalizcd basis. Our matrix 
elements and resulting atomic energies for small cases were validated against an prior, independent code (M. Bromley, 
private communication). 

III. RESULTS 

The following sections give descriptions of the results obtained as well as tables and example figures. The computed 
ground state energies for HF, PHF, and RPA are compared to CI. The ionization potentials and electron affinities 
are compared to experimental values [42j . 

A. Ground State Energy 

The computed ground state energies for the neutral atoms studied are in Table IIVI CI is exact for purposes of 
comparison. As required by variational theory, Hartree-Fock gives an upper bound to the ground state energy and 
projected Hartree-Fock improves upon it. The random phase approximation computes lower energies than CI in all 
of the calculations we performed. There are no CI calculations for several the largest basis sets in combination with 
the largest atoms because the computational requirements exceeded our available resources; keep in mind that our 
'approximate' results were computed in full CI spaces, as we allow arbitrary deformation of the single-particle states. 
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TABLE IV: Ground State Energy (hartree) 



Atom 


space 


HF 


PHF 


RPA 


CI 




CVB1 


-7.4327 


-7.4327 


-7.4720 


-7.4700 


Li 


CVB2 


-7.4327 


-7.4327 


-7.4751 


-7.4725 




CVB3 


-7.4327 


-7.4327 


-7.4774 


-7.4741 




CVB1 


-14.5732 


-14.5758 


-14.7491 


-14.6558 


Be 


CVB2 


-14.5733 


-14.5778 


-14.7271 


-14.6589 




CVB3 


-14.5733 


-14.5780 


-14.7315 


-14.6615 




CVB1 


-24.5239 


-24.5352 


-24.6852 


-24.6391 


B 


CVB2 


-24.5330 


-24.5353 


-24.6965 


-24.6458 




CVB3 


-24.5331 


-24.5355 


-24.7009 


-24.6482 




CVB1 


-37.6700 


-37.7052 


-37.8453 


-37.8202 


C 


CVB2 


-37.6715 


-37.7061 


-37.8661 


-37.8326 




CVB3 


-37.6717 


-37.7072 


-37.8753 


-37.8378 




CVB1 


-54.3443 


-54.4221 


-54.5671 


-54.5517 


K 


CVB2 


-54.3445 


-54.4211 


-54.5971 


-54.5699 




CVB3 


-54.3473 


-54.4236 


-54.6121 


- 




CVB1 


-74.7793 


-74.8267 


-75.0487 


-75.0037 





CVB2 


-74.7822 


-74.8303 


-75.0995 


- 




CVB3 


-74.7829 


-74.8309 


-75.1172 


_ 




CVB1 


-99.4133 


-99.4156 


-99.7283 


-99.6457 


F 


CVB2 


-99.4155 


-99.4188 


-99.7943 






CVB3 


-99.4162 


-99.4195 


-99.8184 






CVB1 


-128.5455 


-128.5455 


-128.9307 


-128.8278 


Ne 


CVB2 


-128.5462 


-128.5463 


-129.0083 






CVB3 


-128.5470 


-128.5470 


-129.0399 





Figure [T] shows the results for carbon in graphical form; the results for all other atoms are qualitatively similar, 
although the error for all approximations increases as we add more electrons. For example, the error in RPA is only 
a few millihartrees for Li, while it is 100 millihartrees for Ne. 

Note there is no difference between the HF and PHF results for Li and for Ne; this is because the HF states for 
neutral Li and Ne are already in states of good total L; projection cannot change this. 

B. Ionization Potential 

Table [V] shows the single ionization potentials for all methods as well as experimental values [42j for all atoms 
studied. Figure O shows the ionization potential for carbon graphically. 

Because the ionization potential (and in the next section, the electron affinity) is the difference between two ground 
state energies, the nice clean trends of the previous section are lost. CI agrees well with experiment, within a 
few millihartree through C and less than 20 millihartrcc difference through Ne. Increasing basis size brings better 
agreement with experiment, though of course we note we have left out relativistic corrections. 

HF consistently underestimates the ionization potentials, with PHF closer for Li through N (with HF and PHF 
identical for Li as the neutral Li and Li I HF states both have good angular momenta) but farther away for O, F, and 
Ne. 

RPA results are the sporadic, sometimes better than HF / PHF and sometimes; again, this is because errors in the 
approximation do not cancel but exacerbate. Increasing the basis size does not always lead to better estimates; they 
get worse for O, F, and Ne. 
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FIG. 1: (Color online.) Carbon ground state energy. 



TABLE V: Ionization Potential (hartree) 



Atom space HF PHF RPA CI 


Exp [42] 


CVBl 
Li CVB2 
CVB3 


0.1963 0.1963 0.1972 0.1969 
0.1963 0.1963 0.1974 0.1971 
0.1963 0.1963 0.1982 0.1978 


0.1982 


CVBl 
Be CVB2 
CVB3 


0.2957 0.2993 0.4324 0.3407 
0.2958 0.3004 0.4072 0.3411 
0.2958 0.3006 0.4090 0.3416 


0.3427 


CVBl 
B CVB2 
CVB3 


0.2955 0.2978 0.1908 0.3002 
0.2955 0.2978 0.1949 0.3030 
0.2955 0.2979 0.1974 0.3041 


0.3051 


CVBl 
C CVB2 
CVB3 


0.3736 0.4063 0.3884 0.4074 
0.3748 0.4068 0.3973 0.4119 
0.3750 0.4074 0.4000 0.4133 


0.4140 


CVBl 
N CVB2 
CVB3 


0.4860 0.5166 0.5218 0.5259 
0.4867 0.5178 0.5312 0.5325 
0.4870 0.5168 0.5343 


0.5383 


CVBl 
O CVB2 
CVB3 


0.4857 0.4320 0.5201 0.4787 
0.4852 0.4335 0.5350 
0.4859 0.4340 0.5408 


0.5007 


CVBl 
F CVB2 
CVB3 


0.6241 0.5674 0.6670 0.6216 
0.6215 0.5652 0.6792 
0.6220 0.5674 0.6858 


0.6405 


CVBl 
Ne CVB2 
CVB3 


0.7258 0.7235 0.8031 0.7769 
0.7220 0.7186 0.8128 
0.7227 0.7192 0.8206 


0.7928 



T 



CI 

■ PHF 

■ — ■ — RPA 

■ — ■ ■ - HF 
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FIG. 2: (Color online.) Carbon ionization potential. 



C. Electron Affinity 



Not all of the second row atoms can bind an additional electron; we only present those cases which do so experimen- 
tally: lithium, boron, carbon, oxygen and fluorine. Table PVTl shows the electron affinities computed by each method 
and the respective experimental values [42j]. Our results are comparable to other work in the field |47j . Figure [3] shows 
the electron affinity of carbon in graphical form. 

The quality of agreement is generally poor, even for CI. In retrospect this is not surprising: we have already seen 
errors for all methods grow rapidly with the number of electrons, and given the smallness of second row electron 
affinities the results are unsurprising if disappointing. Many of the affinities even have the wrong sign. Adding an 
electron to the atom increases the runtimes and memory requirements over the neutral atom calculations, therefore, 
only the smallest basis sets for all but lithium could be used. Without the larger calculations, a clear trend is difficult 
to establish. It appears that the CVB1 basis set is fairly poor when computing electron affinities, sometimes giving 
negative results. However, based on the CVB2 calculations for boron and carbon, the CI electron affinities might 
converge on the experimental values as basis size increases, similar to what is seen with the ground state energy and 
ionization potential calculations. 

All three of our approximations perform poorly, trending neither towards the CI result nor towards experiment. 
HF and PHF are unable to provide consistent results across either atoms or basis sets and generally tend to be too 
far from experimental values to be useful; nor do they show consistent trends as the basis size is increased. RPA 
performs better, but our results do not suggest convergence with basis size. Because it is not clear if this is a fault of 
the approximation or of the basis for CI (which again is our numerically exact benchmark) it would be good, if and 
when possible, to continue this investigation with even larger bases. 



D. Optimization of the space 



We utilized the CVBn STO-basis sets , which had been arrived at through optimization using both Hartree-Fock 
and CI calculations Q. But CI calculations, particularly in large dimension spaces, are tedious and time-consuming. 
Furthermore, a basis optimized for Hartree-Fock calculations may not be optimized for CI calculations. Therefore it 
might be useful to see how well our approximate methods can work as a proxy for full CI when changing the basis. 
In other words, if we vary the basis, does a minimum in CI have a corresponding minimum in one or more of our 
approximate methods? 

We carried out an example optimization by altering all of the parameter values in the STO for carbon CVB1 by a 
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TABLE VI: Electron Affinity (hartree) 



Atom space HF PHF RPA CI 


Exp [421 


CVBl 
Li CVB2 
CVB3 


-0.0152 -0.0152 0.0780 0.0111 
0.0978 0.0050 0.0804 0.0226 
-0.0020 0.0056 0.0768 0.0226 


0.0227 


CVBl 
B CVB2 
CVB3 


-0.0379 0.0163 -0.0302 -0.0141 
-0.0222 -0.0413 0.0009 0.0068 
-0.0205 -0.0004 0.0455 


0.0102 


CVBl 
C CVB2 
CVB3 


-0.0126 0.0081 0.0186 0.0178 
0.0057 0.0203 0.0591 0.0446 
0.0060 0.0192 0.0616 


0.0464 


CVBl 
CVB2 
CVB3 


-0.0110 -0.0564 0.0399 -0.0047 
0.0089 -0.0302 0.1092 
0.0148 -0.0302 0.1148 


0.0537 


CVBl 
F CVB2 
CVB3 


0.0133 0.0090 0.0890 0.0600 
0.0437 0.0405 0.1706 
0.0432 0.0399 0.1748 


0.1250 



0.12 



0.08 - 



0.04 - 



0.00 - 



-0.04 



CI 

— PHF 
■ — RPA 

■ - HF 

■ ■ ' Experiment 



_L 



CVBl 



CVB2 
Basis Set 



CVB3 



FIG. 3: (Color online.) Carbon electron affinity. There are no CI calculation for the CBV3 sets because the 
computational requirements exceeded our available resources. 



scaling factor ranging from 0.5 to 1.5. Table IVHl gives the results of a quadratic fit to each computational method. 
The data show that RPA has a similar minimum to CI. Figure |4] shows these results graphically. 

This optimization suggests that RPA is the most desirable method of those studied for basis set optimization. It is 
much faster than CI and clearly tracks the ground state energy as basis parameters are altered. PHF is a more useful 
upper bound to the ground state energy than HF. but cannot track the CI energy well enough to be used as a proxy 
for CI. 
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TABLE VII: Carbon Minimum Energy Scale Factor 

Method Scale Factor 
CI 0.981 
RPA 0.989 
HF 0.892 
PHF 0.908 




37 g I i I i I i I i I i I i I i I i I i I i I i I i I 

0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 

Scale Factor 

FIG. 4: (Color online.) Carbon optimization. As the STO basis parameters are scaled by a scale facter between 0.5 
and 1.5, the random phase approximation tends to track CI very well. Hartree-Fock and projected Hartree-Fock 
follow each other but not CI. This indicates that RPA can be useful as a proxy for CI during basis parameter 

optimization. 



IV. CONCLUSIONS AND FUTURE WORK 

This paper is an extension of previous work comparing the random phase approximation against CI calculations 
in a CI basis pj, testing how good is an approximation against a numerically exact result. By comparing against a 
numerically exact result, we eliminate any errors due to choice of model space, interaction, etc. As part of this process 
we find eigenvalues of very large CI matrices. To simplify matters we used an uncoupled basis, which allowed us to 
recomputed the many-body matrix elements on the fly in a highly efficient manner. 

Not surprisingly, we found in absolute error RPA gave the best approximation to the ground state energy. RPA is 
not variationally bounded, and in the cases we chose it systematically overestimated the binding energy; this contrasts 
from previous, similar work for nuclei, where RPA both over- and underestimated the binding energy [l|. Because 
the overestimate was not uniform, HF and PHF calculations sometimes gave accidentally better estimates of the 
ionization potentials and electron affinities. 

Perhaps most usefully, we found that for an example fixed atomic system, RPA tracks the CI binding as we varied 
parameters for the basis, suggesting RPA might provide an efficient proxy to speed up optimization of a basis. This 
will be explored further in future work. 

In principle we can also compute and compare transitions and strength functions, as has been done for nuclei @], 
and static and dynamic observables such as static moments, polarizabilities, etc. This we also leave for future work. 

We thank Michael Bromley for many valuable conversations and suggestions, and in particular for aid in validating 
our codes. The U.S. Department of Energy supported this investigation through grants DE-FG02-96ER40985 and 



12 



DE-FC02-07ER41457. 



[1] I. Stetcu and C.W. Johnson, Phys. Rev. C 66, 034301 (2002). 

[2] I. Stetcu and C.W. Johnson, Phys. Rev. C 67, 044315 (2003); I. Stetcu and C.W. Johnson, Phys. Rev. C 69, 024311 
(2004). 

[3] J. T. Staker and C. W. Johnson, to be published. 

[4] I. Lindgren and J. Morrison, Atomic many-body theory, 2nd. ed. (Springer- Verlag, Berlin, 1985). 

[5] P.J. Brassard and P.W.M. Glaudemans, Shell-model applications in nuclear spectroscopy (North-Holland Publishing Com- 
pany, Amsterdam 1977). 

[6] P. Ring and P. Shuck, The nuclear many-body problem, 1st edition, Springer- Verlag, New York 1980. 
[7] I. Ema et al, J. Comput. Chem. 24, 859 (2003). 

[8] D. B. Cook, Handbook of computational quantum chemistry, (Oxford University Press, Oxford, 1998). 
[9] C. Froese Fischer, T. Brage and P. Jonsson, Computational atomic structure: an MCHF approach (Institute of Physics 
Publishing, Bristol, 1997). 

[10] F. Jensen, Introduction to computational chemistry, 2nd ed. (John Wiley & Sons Tld, West Sussex, 2007). 
[11] P.-O. Lowdin, Phys. Rev. 97, 1474 (1955) 

[12] CD. Sherrill and H.F. Schaefer III, Adv. Quant. Chem. 34, 143 (1999). 

[13] B. A. Brown and B. H. Wildenthal, Annu. Rev. Nucl. Part. Sci. 38, 29 (1988). 

[14] I. Shavitt, Mol. Phys. 94, 3 (1998). 

[15] A. Hibbert, Comp. Phys. Comm. 3, 141 (1975); A. Hibbert, Comp. Phys. Comm. 35, C-298 (1984). 

[16] F.A. Parpia, C. Froese Fischer, and I. P. Grant, Comp. Phys. Comm. 94, 249 (1996); F.A. Parpia, C. Froese Fischer, and 

I. P. Grant, Comp. Phys. Comm. 175, 745 (2006). 
[17] W. E. Ormand and C. W. Johnson, to be published. 
[18] E. Caurier and F. Nowacki, Acta Phys. Pol. B 30 (1999) 705 
[19] N. C. Deb and A. Hibbert, At. Data and Nucl. Data Tables, 95, 184 (2009). 
[20] G. Correge and A. Hibbert, At. Data and Nucl. Data Tables, 86, 19 (2004). 
[21] G. P. Gupta and A. Z. Msezane, At. Data and Nucl. Data Tables, 89, 1 (2005). 
[22] M. H. Chen, K. T. Cheng, and W. R. Johnson, Phys. Rev. A 64, 042507 (2001). 
[23] M. W. J. Bromley and J. Mitroy, Phys. Rev. A 81, 052708 (2010). 
[24] J. Olsen, P. j0rgensen, and J. Simons, Chem. Phys.Lett. 169, 463 (1990). 

[25] R. R. Whitehead, A. Watt, B. J. Cole, and I. Morrison, Advances in Nuclear Physics (Plenum, New York, 1977), Vol. 9, 
p. 123. 

[26] E. R. Davidson, Comput. Phys. 7, 519 (1993). 

[27] M. L. Leininger, C. D. Sherrill, W. D. Allen, and H. F. Schaefer III, J. Comp. Chem. 22, 1574 (2001). 
[28] A. W. Weiss, Phys. Rev. 122 (6), 1826-1836 (1961). 
[29] C. F. Bunge, Phys. Rev. 168, 92 (1968) 

[30] F. Saskai and M. Yoshimine, Phys. Rev. A 9 (1), 17-25 (1974). 

[31] C. F. Bunge, Phys. Rev. A 14 , 1965 (1976); At. Data Nucl. Data Tables 18, 126 (1976). 
[32] D. Silver, S. Wilson, and C. F. Bunge, Phys. Rev. A 19 (4), 1375-1382 (1975). 
[33] O. Jitrik and C. F. Bunge, Phys. Rev. A 56 (4), 2614-2623 (1997). 
[34] C. F. Bunge, J. Chem. Phys. 125, 014107 (2006). 

[35] C. X. Almora-Diaz, and C. F. Bunge, Int. J. Quantum Chem. 110 (15), 2982-2988 (2010). 
[36] P.-O. Lowdin, Phys. Rev. 97, 1490 (1955) 
[37] S. Lunell, Phys. Rev. 173, 86 (1968). 

[38] B. Laskowski and S. Lunell, Int. J. Quantum Chem. 9 (S9), 175-182 (1975). 

[39] R. D. Lawson, Theory of the nuclear shell model (Clarendon Press, Oxford, 1980). 

[40] A. R. Edmonds, Angular momentum in quantum mechanics (Princeton University Press, New Jersey, 1996). 

[41] G. H. Lang, C. W. Johnson, S. E. Koonin, and W. E. Ormand, Phys. Rev. C 48, 1518 (1993); E. Y. Loh, Jr. and J. E. 

Gubernatis, in Electronic Phase Transitions, edited by W. Hanke and Yu. V. Kopaev (Elsevier Science Publishers B.V., 

New York, 1992). 

[42] David R. Lide, CRC Handbook of Chemistry and Physics 85th edition, (CRC Press Inc, Florida 2005). 

[43] I. Stetcu, PhD. thesis, Louisiana State University 2003 (unpublished). 

[44] H. Hotop and W. C. Lineberger, J. Phys. Chem. Ref. Data 4 (3), 539-576 (1975). 

[45] H. Jiang and E. Engel, J. Chem. Phys. 127, 184108 (2007). 

[46] A. Graneis et al., J. Chem. Phys. 131, 154115 (2009). 

[47] L. A. Cole and J. F. Perdew, Phys Rev A 25, 1265 (1982). 



